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Abstract 

We present new necessary and sufficient conditions for the existence of 
fixed points in a finite system of coupled phase oscillators on a complete 
graph. We use these conditions to derive bounds on the critical coupling. 

1 Introduction 

The phenomenon of synchronization arises in a wide variety of application areas 
across neuroscience, biology engineering and physics [51 [T71 [21 EE [5] . As such, 
the identification and study of structures and mechanisms that support the 
onset of synchronized behaviour is a key issue in the theory of interconnected 
dynamical systems. In particular, there has been a great deal of interest across 
the mathematics, physics and engineering communities in the development and 
analysis of simple mathematical models of synchronization [^ EE EE US [55] . 

To date, one of the most widely-studied frameworks for the analysis of 
synchronization is the so-called Kuramoto model of phase coupled oscillators 
[TOl [23] . In fact, this model has been used in numerous applications in the 
chemical and biological sciences, and its basic properties have been analysed 
using a combination of numerical and analytical techniques 11, 23, 2HQ]- The 
basic Kuramoto model is comprised of a system of coupled oscillators, which 
may have different natural frequencies, where the coupling between two oscilla- 
tors is given by a weighted sinusoidal function of the difference of their phases. 
The weights used in the model are typically taken to be the same for all pairs of 
oscillators and are given by the ratio of a fixed parameter, the coupling strength, 
to the network size. 

The aspect of the Kuramoto model that has attracted most attention to date 
is the manner in which the onset of synchronization depends on the strength 
of coupling between the oscillators. For instance, at very low values of the 



coupling strength, little or no synchronization is observed. As the coupling 
strength is increased, some partial synchronization appears in the network up 
to a threshold value of the coupling strength, referred to here as the critical 
coupling, at which fully synchronized behaviour emerges [H1I2]- The mechanism 
of (de)synchronization in finite populations of oscillators has been described in 
considerable detail [T31 [T2]- In particular, when the coupling strength drops 
below its critical value and as it continues to decrease, the system undergoes a 
series of so called frequency-splitting bifurcations. At each such bifurcation, the 
ensemble of oscillators subdivides into smaller and smaller groups of oscillators 
with identical average frequency, until eventually all oscillators oscillate at their 
own intrinsic frequency. A detailed analysis of this behaviour for a system with 
three oscillators was given in [13] . While the aforementioned contributions focus 
on the behaviour of the system in the sw&critical coupling regime, the present 
paper studies globally phase-locked solutions, which by definition only exist in 
the supercritical coupling regime. 

Another aspect of the Kuramoto model to have attracted attention recently 
is the emergence of phase chaos [141118) in systems of dimension four and higher. 
A generic feature of coupled oscillator systems, phase chaos in the Kuramoto 
model is most prominent in systems with relatively low dimension (comprising 
between ten and fifteen oscillators) [14] . Again this phenomenon can only exist 
in the subcritical coupling regime, and we shall not further consider it here. 

In the original Kuramoto model, it is assumed that all pairs of oscillators in 
the network are connected with the same coupling strength [10]. This type of 
coupling is referred to as 'all-to-all' coupling and corresponds to a network in 
which the underlying graph is complete [4] . Extensions of the Kuramoto model 
to lattices [8] and rings [20] have also been considered, and more recently the dy- 
namics of coupled oscillators on networks with small- world [S7J [7] and scale- free 
[15] topologies have started to attract a lot of interest. More generally, there 
are many fundamental questions relating to the interplay between a network's 
topology and dynamical processes taking place on it which are still unanswered. 
The work described in [5] , which proposes an extension of group-based symme- 
try, using the so-called groupoid formalism, as a means of classifying possible 
behaviours for networked dynamical systems, is particularly noteworthy in this 
context. 

Many of the recent results concerned with the dynamics and synchronization 
of coupled oscillators have either been based on numerical simulations or else 
have been derived for the limiting case of networks of infinite size. In contrast, 
relatively few rigorous results are available for finite-size networks [53J[3]. In this 
paper, we shall be concerned with synchronization in finite systems of coupled 
oscillators. Specifically: we shall establish (new) necessary and sufficient con- 
ditions for the existence of fixed points in a finite system of coupled oscillators 
(see also [33] [22]); compute bounds on the critical coupling strength for such 
systems; and provide insights into the number of fixed points possible under 
strong coupling. Our analysis is in the spirit of the work presented in J9[ [3] , 
and places particular emphasis on the existence of fixed points. Of course, the 
stability of such fixed points is also a topic of great interest, and has been con- 
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sidered in [HI HH [201 13] • However, we shall not explicitly address the question 
of stability in the current paper. 

The outline of the paper is as follows. In Section 2, we introduce the Ku- 
ramoto model, and review some of its basic properties. Here, we also give a 
formal definition of the critical coupling, which is essentially the lowest value of 
the coupling strength for which fixed points exist. In Section 3, we show that 
fixed points will always exist for sufficiently strong coupling (essentially prov- 
ing that the critical coupling is a finite number), and then, in Section 4 provide 
lower bounds on the critical coupling. Section 5 contains necessary and sufficient 
conditions for the existence of fixed points, which are then used in Section 6 to 
describe an algorithm for computing the critical coupling. Section 7 contains a 
numerical example to illustrate the results of the paper and finally, in Section 8 
we present our concluding remarks. 

2 Mathematical Preliminaries and The Kuramoto 
model 

2.1 Basic Notation 

Throughout the paper, R (C) denotes the field of real (complex) numbers, R w 
(C N ) denotes the vector space of all iV-tuples of real (complex) numbers, and 
R NxN (C NxN ) denotes the space of N x N matrices with entries in R (C). i is 
used to denote the complex number satisfying i 2 = —1. For a vector x € M. N , 
Xi denotes the i th entry of x. Also, Iff denotes the vector in M. N , all of whose 
entries are equal to one. 

We shall use V to denote the projection matrix in ~R NxN given by 

m ■■= w ! J . , , i,j=i,...,N, a) 

and V^R^ shall denote the image of M. N under V. Formally, 

N 

VR N :— {x G : ^^=0}. 

j=i 

2.2 The Basic Kuramoto Model 

The basic Kuramoto model of phase-coupled oscillators under the assumption 
of all-to-all coupling is given by 

k N 

Bi = ^ + -^sin(0, -0i), i = l,...,N. (2) 

Here, 8i(-) € R (S 1 ) and Ui € R respectively denote the phase and intrinsic (or 
natural) frequency of oscillator i, and the constant k € R+ is a global coupling 
coefficient. 
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This model can be described more compactly in vector notation as 



9 = u + kf{6) (3) 

where 9(t) := {Oi{t), ■ ■ ■ , 0jv(i)), w := . . . ,ojn), and the mapping / : R N 
is given by 



f(0 = M))> 

i w 

*® : = ^E sin &-&) !<<<^- ( 4 ) 



The assumption of all-to-all coupling is naturally very restrictive, and ought to 
be relaxed in order for this work to be more directly applicable to the modelling 
of biological systems, or most engineering systems for that matter. Work to- 
wards this end is underway and we hope to be able to present some results in 
the near future. Meanwhile, in this paper, we shall focus exclusively on configu- 
rations with all-to-all coupling. First of all, we recall some fundamental notions 
in the theory of synchronized oscillators. 

2.3 The order parameter 

Let B denote the complex unit disc {z £ C : \z\ < 1}. Then define r : R N D, 
by: 

1 N 

r & ■= ^E e ^- ( 5 ) 

Let r _1 (z) := {£ £ R N : r(£) = z} denote the preimage of r, and note that the 
prcimage is nonempty for all z £ D provided N > 2. We introduce the notation 
TZq := r _1 (0). Then, for £ £ R w , we may express r(£) in polar coordinates: 

y ' [0 t£K K ' 

Here, R : R N [0, 1] and tp : R N \lZo >—>■ [0, 2ir) arc respectively defined as 

R(0 ■- 



-^f>to)j +^Ecos(0)j , (7) 



and 



^(0 := arctan — ^ — . (8) 
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Figure 1: The order parameter r(9) :— R(8)e 1 ^^ is defined as the centroid 
(closed circle) of the set of unit vectors (open circles) associated with the phases 
of the oscillators. 

The following properties of the maps R(-) and ip(-) follow immediately from 
Eqn. ©: 

= \e lc \R(H) = R(t) V£eR"; (9) 

+ cljv) = + c mod 2tt G K^V^o. (10) 

In the physics literature, r(-) is known as the order parameter, and is used to 
characterize the amount of order or synchronization in the system ([2]). The 
idea is to think of the phase 6j of oscillator j as a unit vector e l8j in C; the 
order parameter then corresponds to the geometric centroid of the set of vectors 
{e i9j : j = 1, . .., N}, as illustrated in Figure [1] The magnitude of the order 
parameter, given by R(9), serves as a measure of the order in the system, in the 
sense that the closer the vectors are to being perfectly aligned, the closer R(9) 
is to its maximal value 1, while vectors that are far from alignment will give rise 
to values of R(0) significantly smaller than 1. 
It follows from © that for £ e R N \K , 

m = e -^V(0 (11) 

= ^E^"" (4)) - (12) 

Equating real and imaginary parts in (|11| . we immediately see that for £ € 
R N \K : 

1 N 

N 

5>in(^)-&)=0. (14) 



i?(e + cijv) 



i " 

TV E 
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Both of these identities shall prove useful throughout the paper. 

Before proceeding, note that the function / : M. N h— >■ M. N given by ([4]) can be 
written in terms of the functions R(-) and ip(-) as 

f (£) ■= /^(OsinWO-fi) CGR^TZo , , 

for 1 < i < AT. 



2.4 Fixed points and global phase-locking 

Let (ui) denote the sample mean of the natural frequencies, (ui) := Ylj=i w j- 
Similarly, let (9(t)) denote the mean phase of a solution of © at time t. In 
general, (lo) and will be non-zero. However, we shall now show that 

for the study of phase-locked solutions of we may assume without loss of 
generality that (u>) — 0, (0(f)) = for t > to- This helps to simplify the analysis 
of phase- locked solutions of (J2j) , as it allows us to transform the problem into a 
question of fixed point existence for a lower-dimensional system. 
Consider the new coordinates 



Xi(t) := 9i(t)-(d(t)), i = l,...,N. (16) 

Then x(t) := V6(t). Similarly, define := Vu>. In the new coordinates, the 
system dynamics are given by: 

x = n + kf(x), x(t)eVR N , (17) 

where /(•) is defined in ((4]). 

The key point here is that as (Q) = 0, VM. N is invariant under (fl~7|) . To 
avoid confusion, we shall use x(t) to denote solutions to the system (Ti~T|) on 
VR N , while 9{t) shall be used to denote solutions to the original Kuramoto 
system ^ on M. N . Our main concern for the remainder of the paper is to find 
conditions on k and f2 under which the system (|17p has one or more fixed points 
in the sense of the following definition. 

Definition 1 (fixed point) Given u> £ M. , let f2 := Vu>. We say that x S 
VR N is a fixed point ( of the system Jj7| ) ) if 

kf{x) = -n. (18) 

There is a natural correspondence between fixed points of (|17[) and phase- 
locked solutions of ([2]). In fact, for every fixed point x* € VR N there is a 
1-dimensional manifold M := {0 e R N : 9 = x* + (co)t, t e R} that is invariant 
under the original system dynamics ([2|). More precisely, let x* be a fixed point 
and let 9° e R N be such that V9° = x* . Then the solution 9(t) of the system flS) 
with initial condition 6 (to) = 9° satisfies 

i (t)-0 j (t)=6P i -d° j (19) 
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for all t > to and all In other words, a fixed point in the sense of Dcfini- 

tion[T]corresponds to a situation in which each oscillator is phase-locked to every 
other and moves at constant speed Oi = (u>) . We shall refer to this phenomenon 
as global phase-locking. In the literature, it is also known as full (or complete) 
synchronization. See also [§]■ 

2.5 Critical coupling 

We next define the notion of critical coupling , which is central to the work of 
the rest of the paper. Essentially, the critical coupling is the smallest k for which 
the system (|T71) has at least one fixed point. Formally, we have the following 
definition. 

Definition 2 Given uj <E M N , let il := Voj. We define the critical coupling, k c , 
as follows: 

k c := inf {k e M+ : 3x e VR N s.t. kf(x) = -0} . (20) 

k 

Note that this definition of the critical coupling, which is equivalent to that 
of Kl in [9] , does not coincide with the traditional notion used in the physics 
literature. Indeed, the traditional notion of critical coupling is defined in terms 
of the lowest value of k for which there exists at least one solution x(t), t > 
to, and a constant c £ (0,1] such that R(x(t)) = c for all t > to (so called 
stationary or steady solutions [23]). Note that these solutions are not necessarily 
fixed points, although, in finite dimensions, the probability of finding stationary 
solutions that are not fixed points is vanishingly small. In his original analysis, 
Kuramoto showed that in the limiting case when N tends to infinity, stationary 
solutions always exist for large enough k, provided the distribution of natural 
frequencies is symmetric. Our definition, although more restrictive in a sense, 
does not impose any restriction on the shape of the distribution of natural 
frequencies other than that it should have compact support. In fact, it follows 
from the result of Lcmma[4]below that, if the distribution of natural frequencies 
does not have compact support, then the critical coupling will exceed any finite 
number with probability tending to 1 as N tends to infinity. In this paper we 
shall therefore focus on distributions with compact support. 

3 Existence of fixed points under strong cou- 
pling 

In this section we shall show that, provided the distribution of intrinsic frequen- 
cies has compact support, the critical coupling given in Definition [2] is always 
finite. Following on from this, in the next section, we shall derive a number of 
lower bounds for the value of the critical coupling. There are two steps in the 
derivation given here: first, we characterize the fixed points of the homogeneous 
system 

x = kf(x). (21) 
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Then, in the second step, we use a perturbation argument to show that for every 
fixed point of the homogeneous system, we can find an open set containing it, 
such that, under strong enough coupling, the original system (fTT|) has a unique 
fixed point on this set. As a first step, the following lemma characterises the 
fixed points of the homogeneous system. 

Lemma 1 Let /(•) be given by |7J) and £ £ K . We have that /(£) = if and 
only if one or both of the following conditions are satisfied 

(a) R(£) = 0; 

(b) Bin(& - £j) = for all (i,j). 

Proof: Sufficiency of conditions (a) and (b) follows from (fTS)) and (|U) re- 
spectively. To prove necessity, suppose /(£) = and ^ (if R(£) = we 
are done). It follows that sin(-0(£) — ^) = for all i. This implies that there 
exist integers fcj € Z, i = 1, 2, . . .N such that "0(C) ~ & — ^i 71 " f° r au an d we 
have that £j — £j = (fcj — fci)7r. We conclude that sin(^ — £j) = for all (i,j). 

m 

Remark 1 It is not hard to see that conditions (a) and (b) in Lemma\]\ are mu- 
tually exclusive if and only if the dimension N is odd. We shall prove necessity. 
Suppose conditions (a) and (b) both hold and suppose furthermore that N is odd. 
Then for all (i,j) we have that either cos(£i — = 1 or cos(f; — £j) = — f . 

We write R 2 (() = fa ^ . cos(& -Q = fa(N + 2 V ; C os(& - 0)) - «nce 

— by assumption, it follows that 2^^ j>i cos (£i — £j) = —-/V. TTie /e/t 
/land side evaluates to an even integer. By assumption, the number on the right 
hand side is odd. We arrive at a contradiction and conclude that if N is odd, 
conditions (a) and (b) cannot both hold. 

Next we shall prove that the fixed points of our TV-dimensional system (| 1 T|) 
can be found by solving a system of N — 1 equations in as many variables. We 
have the following result: 

Lemma 2 Let p E {1, . . . , N} and let x* G VR N . Then x* is a fixed point 
of JiTf ) if and only if kfi(x*) — — f2, for i ^ p. 

Proof: The proof of necessity is trivial. To prove sufficiency, recall that 

JV 

^2(n j +kf j (x))=0 for all x S R N . (22) 

i=i 

Now suppose kfi(x) = — for all i p. Then it follows from that 
+ kf p (x) = 0. In other words, it follows that kfi(x) = —Hi for all i. We 
conclude that a; is a fixed point. ■ 
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Let x* £ VM. N be a fixed point of the homogeneous system (f2"Tj) such that 
R(x*) 0. We shall now show that locally, in a neighborhood of x* , the system 
of equations 



-Oj = kfi(xi,...,Xff-i,-'£,-iXj) 



(23) 



-Ojv-i = kfN-i(xi, ■ ■ ■ ,xn-i, - Ylf=i x j) 



has a unique solution, provided k is large enough. It follows directly from 
Lemma [2] that every solution of (|23|) defines a fixed point and, conversely, that 
every fixed point satisfies (|23|) . We proceed as follows. Let x* £ VTR.^. We 
define the Jacobian J(x*), as follows: 



[Jij{x )] := 



dxj 



(24) 



where i, j = 1, . . . , N — 1. We have the following result: 



Lemma 3 Let /(•) be given by ^ and suppose that x* £ VM. N satisfies f(x) = 
and R(x) ^ 0. Then det(J(x*)) ^ 0. 

Proof: Let x* be a fixed point of the homogeneous system and suppose 
R(x*) ^ 0. Then by Lemma[TJ we have that sin(x* — x*) = for all (i,j), and 
it follows that 

cos(x*-x*) = cosdxj - x* s ) - (x* - X*)) 

= cos(x* — x*) cos(x* — x*) (25) 

for all s and all The claim is that J(x*) is nonsingular. To prove this, we 

proceed as follows. From the definition, we have that 

j , n = f-E^W cos «-<)- 2cos K-^) (26) 

\cos(^-x*)-cos(x^-x*) i?j 
Using the aforementioned identity, setting s = N, we rewrite (j2l)]) . as follows, 

[ (cos(x* - x^) - 1) cos(x* - x^) I 7^ J 

Inspection shows that the rank of J(x*) is invariant under permutations of the 
components of x*. Hence we can assume, without loss of generality, that there 
exists p £ {0, . . . , N — 1}, such that 

cos(x*-x^) = l ^ l ^P (28) 

V 3 N> 1+1, p+l<i<N V ' 
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Under this assumption J (a;*) takes the form 



J(x*) 



A 
C B 



) 



(29) 



where A and i? are square matrices of dimension px p and (N— 1 — p) x (-/V — 1 — p) 
respectively. It follows that J(x*) is nonsingular if and only if A and i? are 
nonsmg ular. Inspection shows that A = (N - 2p)I + 211 T and B = (2p- N)I. 
It follows that A or B is singular if and only if iV = 2p. In case iV is odd, this 
condition is never satisfied. In case N is even this condition, combined with (|25p 
and the fact that R 2 (x) = cos(xi — Xj), implies that R(x*) = 0, which 

contradicts our starting assumption. We conclude that, under the hypotheses 
of the lemma, J is nonsingular. This concludes the proof. ■ 



and note that II has an inverse II 1 that is defined everywhere in VR^. We are 
now ready to state the main result: 

Theorem 1 Let /(•) be given by and x* g VM. N be such that /(#*) = 
and R(x*) ^ 0. Also, let ft £ UR^. Then there exists K G M, and an open 
set U G R^" 1 such that (a) II -1 (x*) is an interior point of U; and (b) for all 
k > K , the system of equations i2St\) has a unique solution on U . 

Proof: Define y* :— II -1 (x*), and let g : IR^ -1 i— > l^" 1 be given as 



9i(y) ■= fi(Vl, ■ ■ ■ ^JV-i^Ejti 1 Vj), i = h-..,N-l. Note that g(y*) = 0. 



Also, by Lemma |31 we have that det(^(j/*))| ^ 0. Under these conditions, 

the Inverse Function Theorem says that there exists an open set U C R^ -1 
containing y* such that g\ v : U i— > g(U) is a diffeomorphism. By continuity 
(and bijectivity) of g^ 1 there exists S > such that for all z € R w_1 satisfying 
||z| < 5, the equation g(y) — z has a unique solution on U. Now let z be given 
as Zi := —VLi/k. Since, by assumption, max^ \Q,j\ < oo, it follows that, provided 
k is large enough, the system of equations {kgi{y) — — Q : i = 1, . . . ,N — 1} has 
a unique solution on U. This concludes the proof. ■ 
As alluded to earlier, there is a unique correspondence between solutions 
of (|23]) and the fixed points of the system (fTT)) . Indeed, by Lemma [2] we have 
that if y is a solution of ([23]) . then n~ 1 (y) is a fixed point, and conversely, if y is 
a fixed point, then H(y) is a solution of (|23|) . Thus, an immediate consequence 
of Theorem Q] is that for large enough k, the system (f23| will have at least one 
fixed point. In other words, Theorem [T] tells us that the critical coupling, fc c , is 
always finite. 

Note furthermore that the proof of Theorem Q] does not require detailed 
knowledge of the coupling function g and that, as such, its applicability is not 
restricted to networks with all-to-all coupling. To illustrate this, consider the 



Let IT : R 



I— > VM. , be given as 




fori = l,2,...JV-l; 
for i = N. 



(30) 
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case of a 4-cycle, where g is given as 



9i{y) = ^sin(?/2 - Vi) + ^sin(-2yi - y 2 - y 3 ) 
92(y) = -rsm(y 3 -y 2 ) + -sm(y 1 -y 2 ) 

33(2/) = j sin(-2y 3 - yi - j/a) + - sin(y 2 - 2/3)- (31) 

We have that ff (0) = and det(||(0)) =-1^0. This implies that for k 
large enough, the system of equations {kgi(y) — ~Qi : i = 1, 2, 3} has a unique 
solution on some open set containing the origin. 

Lastly, note that continuity of g^ 1 implies that the fixed points of the original 
system (|T?|> will converge to the fixed points of the homogeneous system (|2Tj) 
as k tends to infinity. 



4 Lower bounds on the critical coupling 

In the previous section we showed that the critical coupling is finite, provided 
the oscillator's intrinsic frequencies are finite. In the present section we shall 
investigate in more detail the relation between the distribution of intrinsic fre- 
quencies and the critical coupling. In particular, we shall derive various lower 
bounds and discuss some of these bound's implications for the system's dynamic 
behaviour. 

First, let us observe that k c (Definition [2]) is lower bounded by the norm 
of fi: 

k c > Halloo := max|a>j — (oS)\. (32) 

i 

This follows trivially from inspection of Eqn. (|17p . In order to derive another 
lower bound, we shall need the following result: 

Lemma 4 Let /(•) be given by Then: 

1. For allxe R N , 

\\f(x)h < v / ^W r ^»; (33) 

2. If N is even, then for every c £ [0, 1] there exists x £ VM. N such that 
R(x) = c and \\f(x)\\ 2 = y/NR?{x) (l - R 2 (x)); 

3. If N is odd, then inequality h33\) is strict for all x £ R N such that < 
R(x) < 1. 

Proof: Part 1. Observe that inequality (|33|) is trivially satisfied when 
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x € 72-0- Suppose therefore that x € K \TZo- Then by definition 



JY 



l/Wlli 



(34) 



where VK 2 -) an d R(x) are the phase and magnitude of the order parameter, pre- 
viously defined in ([5]) and ([7]) respectively. Introducing the shorthand notation 
Zi(x) := cos(ip(x) — Xi), and using (TIB"]) we now rewrite (fM]) . as follows: 



' jv \ 2 JV 

EC 1 

v J =1 / '' =1 



(35) 



To derive the desired inequality we pick a c € [0, 1] and maximize ||/(x)||2 over 
the set {x £ M. N : R(x) = c}. We shall not solve this optimization problem 
directly, but take an indirect route by considering another, easier optimization 
problem, whose solution will then give us an upper bound on the solution to 
the first problem. Then we shall show that, under certain conditions, the two 
solutions coincide. 

To this end, let c S (0, 1] and consider the constrained optimization problem 



OPT 1: 



maximize 



Eli (i-*(*) a 



subject to i J2f=i 



Zj(X 



x £ 



Note that the constraint is feasible for all values of c in the specified interval. We 
shall denote the solution to OPT 1 as si(c). Next consider a second optimization 
problem, 



OPT 2: 



maximize 



subject to ^Y,j=iVj = c j 



y e 



pN 



and let the solution to this problem be denoted as S2(c). We then have that 
s 2(c) > si(c) for all c £ (0,1]. In other words, the solution to OPT 1 is 
upper bounded by the solution to OPT 2. The solution to OPT 2 can be 
found by means of standard Lagrange multiplier techniques. The optimum 
52(c) = N (l — c 2 ), is attained when yi = c for all i. We conclude that 



max II /(a;) II 2 , < Nc 2 (l - c 2 ) , 



and hence, 



|/(x)|| 2 < VNR(x)y/l-R?(x). 



(36) 



(37) 
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for all x e R N . 

Part 2. To prove the second part of the theorem, let c E (0, 1] and note that 
si(c) = S2(c) if and only there exists x € R JV \7?-o such that 



cos(V>(x) — Xi) = c 
for all i. Suppose N is even and let x be given as 

J arccos(c) i = 1, 
1 — arccos(c) * = y 



jv 

2 



, AT. 



(38) 



(39) 



Then SjLi 35 j = 0; an dj by definition, x £ VI 
cos(i/j(x) — Xi) = c for all i. This completes the second part. 

Part 3. To prove the third part, let N be odd and suppose there exists 
x £ M. N such that Condition l|38p is satisfied. Then it follows from the identity 

N 



pN 



Moreover, ip(x) 



0, and 



sin 2 ('0(x) — Xi) + cos 2 (ip(x) — Xi) 
such that sin(^)(x) 



1 that there must exist a e { — 1,1} 



) = diVi — c 2 for all i. By Identity (fT4| . we have that 



J2j sin(t/j(x) — Xj) = 0, which, assuming c ^ 1, implies that J2jLi a j = 0- Brrt 
this cannot be true unless N is even. Thus we arrive at a contradiction and we 
conclude that if N is odd then s 2 (c) > Si(e) for all c such that < c < 1. This 
concludes the proof. ■ 




0.2 0.4 0.6 0.8 1 

R(x) _ 

(a) TV = 3 




0.2 0.4 0.6 0.8 

R{x) 

(b) N = 4 



Figure 2: Scatter plot of -^=||/(x)|| 2 for N = 3 (left panel) and iV = 4 (right 
panel). The phases x were drawn from a uniform distribution. The solid black 
line in both panels is the upper bound R(x)y/1 — R 2 (x). 



Figure [2] illustrates the result of Lemma [4] When N = 4 (even), the lower 
bound is attained at every value of R(x), which shows that the given bound is 
the tightest possible. However, as illustrated in the left panel, when N = 3, the 
bound is never attained except on the set {x € K w : R(x) £ {0, 1}}. It can be 
shown however that in the limit of large N the given bound is arbitrarily tight, 
even for odd N, in the sense that for every c G [0, 1], 

||/(x)|| 2 - v/(2m + l)c 2 (l-c 2 ) 



1 

mm 



»+ 1 :Rfa ! 1=fil */9.m. 4- 1 
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tends to zero as m tends to infinity. 

Lemma [4] has some interesting implications. For instance, it can provide 
insight into the rate at which solutions of a homogeneous system of Kuramoto 
oscillators (© with u>.- L — for 1 < i < N) on R N converge to fixed points, . To 
see this, consider the homogeneous system 

Z : kI T ■ «-> 

where 9q € M. N . We shall compute the time derivative of the magnitude squared 
of the order parameter, L(-) := R 2 (-), and show that this derivative is: (i) non- 
negative along solutions of (|4"0f : (ii) bounded from above by a certain function 
D{t) for every t. We proceed as follows [5]. By definition, 



Using the identity 



it follows that 



dL(9(t)) L(0). L (0), f(f)M , 



89 N 

mmi (4i) 



dL(0(t)) _ 2k nj . ilii: , t 2 



dt N' 

which shows that the time-derivative is positive everywhere, except at the equi- 
libria, where it is zero. It follows that the magnitude of the order parameter 
is a nondecreasing function of time. Based on the observation that the time- 
derivative of L is positive almost everywhere (the set of equilibria having mea- 
sure zero), we formulate the following conjecture, OE]: 

Conjecture 1 For almost all initial conditions 8q, the solution 8(t) to the ho- 
mogeneous system J^0| j has the property that Hindoo R(9(t)) = 1. 

In agreement with Conjecture [TJ one can prove that, for the homogeneous 
system, the global phase- locking manifold M. := {6 E M. N : &i — Oj for all i,j} 
is (locally) asymptotically stable. However, the existence of other invariant 
manifolds, not contained in yVf, implies that M. is not globally asymptotically 
stable. We conjecture that M. is 'almost globally asymptotically stable', in the 
sense that its region of attraction is the entire space minus a set of measure 
zero. 

For our next result, we shall need the concept of a dominating function, 
which is defined as follows: 

Definition 3 Let h, g : R i— ► M. and let X C R be some interval. We say that h 
dominates g onX ifh(t) > g(t) for all t £ X. In that case we call h a dominating 
function for g on X. 

Our next result states that L(9(t)) is dominated by a certain scalar function 
D{t) that depends only on 9q. In order to prove this result, we need the following 
two lemmas. 
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Lemma 5 Let #(•) be a solution of the homogeneous system and suppose 
L{9{t')) = 2kL{9{t')) (1 - L{6{t'))) for some t' G R. Then 

L(6(t)) = 2kL(9{t)) (1 - L{9(t))) for all t > t' . 

Proof: Recall that L(9(t)) = ^\\f{9{t))\\ 2 2 . It follows from the proof of 
Lemma H that ||/(0 = NL(6{t')) (1 - L(9(t'))) for some t' G R if and only 
if one or two of the following conditions hold: (a) L(9(t')) = 0; (b) N is even 
and there exists a permutation 9(t') of 9(t') such that 

N 

cos(0i(i') - = 1, i = l,2,...,- (42) 

AT 

cos(^(t')-^Jv(^)) - 1. i=y + l,...,iV (43) 

If L{9{t')) = 0, we have that = for all t > t' and the result follows 

trivially. Now suppose Conditions (|4*2"j) and (|4*3")) hold. It follows that 9i(t') = 
9j(t') for i, j < y and i, j > y, and hence 



A (cos(^(r)-^(r)) 



iV AT 

o, hi < y ; m > y 



In other words, if (?(•) satisfies Conditions ([42]) and (|43f for some i', it satis- 
fies (J42J) and ((43j) for all i > t'. This concludes the proof. ■ 

Lemma 6 Let h, g : R i— > R fee suc/i i/iat /or every xq G R, i/ie systems 

x — h(x) ( x = g(x) 

x(0) — xq ' \ a;(0) = 

Ziaue unique solutions in C^O, oo). Let £/iese solution be denoted Xh(t;xo) and 
x g (t;xo), respectively. Suppose there exists a, b G R, a^b, such that 

h(x) > g(x) > 

for all x G (a, 6) C R. Then for every x$ G (a, 6), we have that 

x h (t;x Q ) > x g (t;x ) (44) 

for all t £ I , where I C R is defined as 



[0, min t {x h (t; x ) = b}) when {x h (t; x ) = 6} ^ 
[0, oo) otherwise. 



(45) 



Proof: Under the hypotheses of the lemma we have that h(xo) > g(xo), 
and it follows that for small enough t, Xh(t) > x g (t) (omitting the argument Xq 
for notational convenience). Note also that Xh and x g are increasing whenever 
Xh{t\Xo) < b and x g (t;xo) < b respectively. Now suppose there exists t 2 > 
such that x g (t 2 ) > Xhfa) and x g (t 2 ) < b. Then by continuity there exists 
ti < t 2 such that a < Xh(t\) = x g (t\) < b and h{xh{t\)) < g(x g (ti)). Define 
x' := Xh{t\) = x g {ti). It follows that g{x') > h(x'), which contradicts our 
starting hypothesis. We conclude that Xh(t) > x g {t) for all t G I. ■ 
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We have the following result: 

Corollary 1 Let #(•) be a solution to the homogeneous system with initial 
condition 9 (to) = %■ Then 

D(t) := r (46) 

l- e - 2 *(-o)(W-) 

is a dominating function for L (9 (t)) on [to,oo). 

Proof: By Lemma H we have that L(9(t)) < 2kL(9(t)) (1 - L(9(t))) for 
all t. We claim that, on [to, oo), L(9(t)) is dominated by the solution y(t) of the 
ODE 

V = 2fcy(l-y) 



y(to) = L(9 ) (47) 



which is given as 



To prove this, suppose L(9(t)) = 2kL(6(t)) (1 - L(9(t))) for some t > t and let 
denote the smallest such t (in case L(9(t)) < 2kL(6(t)) (1 - L(6(t))) for all 
t > to, the result follows immediately from Lemma [6]). Then by Lemma we 
have that L(9(t)) = 2kL(9(t)) (1 - L(6(t))) for all t > t' , and it follows that 

^(*))lw))=a= 1 _ e _ 2fe( L) (a^y *^ ( 49 ) 

One can easily verify that L(9(t))\ L ^ g ^ t ,^ =a is nondecreasing as a function of a 
for all t > t' (and a € [0, 1]). Now let 2 be an upper bound for L(9(t')). It follows 
that - — _ 2 fc(t-t') / n 'i i s a dominating function for L(9(t)) on the interval [t', oo). 



-7.klt.~t.' \ ^ l-l j 



To compute an upper bound for L(9(t'), we proceed as follows. By definition 
of t', we have that L(6(t)) < 2kL(9(t)) (1 - L(9(t))) for all t < t'. It follows 
from Lemma [5] that 

mt)) * j ,-^-i r ^)-iv ^*<* (50) 

By continuity, we have that 

L(9(t')) < - 



Using the upper bound — 1 , L(e , _ lN for L(9(t')), it follows from (|49[) 

1 _ e -2fc(t -t )^ ^ j 

that 

- w TTXrel IT ' *^ (51) 
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Combining (jST))) and (f5Tj) , we arrive at the desired result. This concludes the 
proof. 

■ 

Figure [3] shows the graph of L(9{t)) and that of the dominating function 
D(t) — Eqn. (|46|) for a particular realization of the initial condition 9q. In this 
example, N — 100 and k = 2. We observe that, in agreement with ConjectureQ] 
the solution converges to a globally phase-locked state, that is L((9(t)) — > 1. 
Note that convergence can be very slow depending on the choice of initial con- 
dition. Indeed, for any Tel and any e > 0, we can find 5 > such that if 
L(6 ) < 5 then L(6(t)) < e for all t < t + T. The upshot of this is that if 
the initial condition 6o is selected by drawing from a uniform distribution and 
the number of oscillators is large, then L(9q) is likely to be small, and as a 
consequence convergence to the stable equilibrium is likely to be slow. In the 
limit case when N tends to infinity, we have that L(9q) tends to zero with prob- 
ability 1 and the time required for L(6(t)) to exceed some given finite threshold 
diverges to infinity. 




1 2 3 4 5 

t -> 



Figure 3: Numerical simulation of the homogeneous system ([40|) with N = 100 
oscillators and coupling coefficient k = 2: time evolution of L(0(t)) := R 2 (6(t)) 
(solid line) and the dominating function D(t) — Eqn. pp]) (dashed line). 



Let a u := y^X)j=i ( w j ~ ( w )) 2 denote the (sample) standard deviation 
associated with the vector of intrinsic frequencies uj. Using Lemma [4] we derive 
another lower bound on the critical coupling, as follows: 

Corollary 2 The critical coupling k c satisfies 

k c > 20-uj (52) 

Proof: Let x* E VM. N be a fixed point of the system (fT7| . Then by definition 
fc||/( a; *)l|2 = | V^oj 1 1 2 = \fNa u and by LemmaHwe have that 

||/(z*)|| 2 < y/Ny/R*{x*)(l-R*{x*)). (53) 
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Figure 4: Graph associated with inequality ([55)1 . For a given value of the ratio 
(fu/fc), the magnitude of the order parameter evaluated at a fixed point 
a;*, must lie within the striped region. 



It is not hard to see that the right hand side of ([53]) is upper bounded by \^/~N ■ 
It follows that 

^ » (54) 



k > 



2a„ 



This completes the proof. ■ 
Note that Corollary[2]is in agreement with the intuition that greater variation in 
intrinsic frequencies requires stronger coupling to achieve global phase-locking. 

Using Lemma [4] we can compute bounds on the value of the order parameter 
evaluated at the fixed points of the system, should they exist. Indeed, suppose 
k > k c , then for any fixed point x* £ VM. N we have that 



y/B?(x*)(l-R 2 (x*)) > 



(55) 



Solving for R(x*) gives 



1 1 



1 1 



(56) 



The graph associated with inequality ([5fi|) is shown in Figure 2] 



5 Necessary and Sufficient conditions 

In the last section, we derived lower bounds for the critical coupling of the system 
(|17p . which provide necessary conditions for the existence of fixed points. We 
next derive conditions that are both necessary and sufficient for fixed points to 
exist, and we shall use these results to describe an algorithm for computing the 
critical coupling in Section 6. 
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Throughout this section, the function / : R N i— > RR N is given by ([U and 
the set T{k, f2) is defined as 

T{k, 0) := {a; € VTR^ : fc/fx) = -0} , fc > 0, fie VR N 

On .F(fc, f2) we introduce a notion of equivalence, as follows. 

Definition 4 (Equivalence on F{k, Q)) Given Q £ VR N and k > 0, let 

x,x' € F(k 7 Q). We say that x and x' are equivalent (x ~ x' ) if R(x) — R(x'). 

To motivate this definition consider the following fact. Let fc, fi be given, and 
let s e Z N be such that YljLi s j = 0- If £ is a fixed point of the system (jTTJ) , 
then x' := x + 2sir is also a fixed point of the system (|17[) and, in addition, 
R(x') = R(x). 

The following theorem provides a necessary and sufficient condition for the 
system (|17[) to have a fixed point, given a particular coupling strength fc, and a 
particular realization of intrinsic frequencies, fi. 

Theorem 2 Let fc > and Q e VR N . Then T(k,Vt) ^ iff there exists 
P € [il|fi||oo, l]cl and a g {-1, 1}^ suc/i t/wif 



3=1 



Moreover, suppose (a ,f3 ) and (a 2 ,/3 2 ) both satisfy |57| ) and let x 1 ^ 2 € F(k,Q) 
be such that 

kPMMx')-^) = "fy ie / 12 > 7 -12 JV (58) 
a,cos(^(x 1 )-^-) > ' je ^>^> 1,^, ...JV ^Sj 



27»en x 1 ~ x 2 iff /3 1 = /3 2 and ^ =i (<H - aj) V 1 - ( ]j* J = 



a 4 x 2 



Proof: Suppose £1 7^ (the case f2 = is easy). Let x* £ VR N be a 
fixed point of (fTT)) . By definition, kf(x*) = —CI, and since fi 7^ 0, we have that 
f(x*) ^ 0, and consequently R(x*) 7^ 0. It follows that 

sm(iP(x*)-x*) = - I ^, i = l,2,...,N. (59) 

Let [3 := R(x*). By dS5J) we have that /3 > |||fi||oo- Recall that for all x € 
M. N \lZo, R(x) can be written as 



3=1 

and let ai be given as 

fl . f-1 ifco S (V>(**)-<)<0; 
I +1 otherwise. 
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Combining (|5^|) , (|6"0)) . and (j6"Tj) , we arrive at 

1 N 

"4i 



A-/? 



(62) 



This proves necessity. To prove sufficiency, let a G { — 1,1}^ be given, and 
suppose [3 > r||^||oo > (again, the case SI = is easy). Then for every c e E, 
the system 



kj3sm(-y l - c) = -fij 
di cos(-y i - c) > 0, 



1,2, 



has a unique solution y* e [— n,ir) N . We pick c such that Xw=i 2/j = 0- Since 



(63) 



Sj-Li sm (?/j + c ) = 0, it follows that 
R(y*) = R(y* + cl) 
From (l6"3l, we have that 



cos(y* + c) 

3=1 



cos(y* + c) = a;y 1 - 
Combining (164)) and (|65p . we arrive at 

N 



k3 



i = l,. 



, AT. 



N ^ 



3=1 



fc/3 



(64) 



(65) 



(66) 



The second part of the theorem follows easily after noting that if x % satisfies (|58l) 
then R(x l ) = [3\ i = 1,2. ■ 
Theorem [2] gives us a necessary and sufficient condition for the equation 
kf(x) = — f2 to have at least one solution for a given value of k. It is not clear, 
however, that there exists a k for which this condition is satisfied. The following 
Corollary provides an easy sufficient condition. 



Corollary 3 Let k > and Q e VR N . Suppose 
1 1 N I 

3=1 ' 



(67) 



for some a € {-1, 1} N . Then T(k, 9) ^ 0. 

Proof: Suppose fl ^ (again, the case O = is easy). Let a £ { — 1, 1}^. 
Define m : [iHfi^,!] i-> R, m (/3) := and n : [£||fi|U, 1] X {-1,1}^ ^ M, 



3=1 



fc/3 



(68) 
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Since f2 ^ we have that m(l) > n(l, a). Now suppose condition ([S7]) is satis- 
fied. Then we have that r7x(-^ || 0|| oo) < Ji(^||f2||oo( a), and by the Intermediate 
Value Theorem there must exist (3* G [^ll^Hoo)!] sucn that m(j3*) = n(f3*,a). 
It follows from Theorem [5] that the system (|T71) has a fixed point. ■ 

Corollary 4 Let il G VR^. Then: (i) the critical coupling k c is finite; (ii) for 
large enough coupling, the system |J7| ) has at least 2 W ~ 1 fixed points. 

Proof: Note that the right hand side of (|67p does not depend on k. Hence, 
it follows that, provided 



j=1 V Ml 

Condition ([67| is always satisfied for large enough k. Furthermore, it follows 
easily that if (|69|) is not satisfied for given a, then it is satisfied for a' := —a. 
This implies (i) that the critical coupling k c is always finite, and (ii) that the set 
A + := {a G { — 1, 1}^ : Eqn. ([M]) is satisfied} contains precisely 2 JV_1 elements 
(counting multiplicity), each of which defines a unique (up to equivalence in the 
sense of Definition [4]) fixed point. This concludes the proof. ■ 

Corollary 5 Let k > and Q G VR N . Then T(k, Q) ^ if and only if there 
exist (3 G [i||fi||<x» 1] suc h that 




k/3 



Proof: The proof of Corollary [3] suggests that if the fixed point equa- 
tion (|57p does not have a solution, then necessarily 





for all (3 G [rllfiHooj 1] and for all a G { — 1, 1}^. Since we have that 
1 N 

-Y 

3=1 

for all a G {— 1, 1} , it follows that the given condition is necessary and sufficient 
for the system (fTTj) to have at least one fixed point. This concludes the proof. 

■ 

The next and final corollary gives us an upper bound on the critical coupling. 
Corollary 6 The critical coupling, k c , satisfies: 

k c < Mo ° (70) 



Proof: Follows directly from Corollary [3J 
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6 An algorithm for computing k c 



In this section we present a bisection algorithm that will allow us to numerically 
evaluate the critical coupling with arbitrary precision. Throughout, we shall 
assume that £1^0. Define X := (HQH^, oo), and let pi : X i— > (0,1], and 
P : X t— * (0, 1] be given as 

Also, define /i(u; k) : X x R + M + , 

Jfe) := i«. (72) 

From Corollary [5] it follows that the critical coupling is the smallest k for which 
the equation P(u) = h(u; k) has at least one solution on X. We have the following 
result. 

Theorem 3 For all Q 6 VR N , Q^O, the equation 




has a unique solution u* G X, and we have that 




Proof: Observe that, by strict concavity of P and linearity of h( ■ ; k), the 
equation P(u) = h(u; k) can have at most two solutions on X for any k > 0. We 
shall now show that, when k = k c , it can have no more than one solution. Since, 
by definition of critical coupling, P(u) — h{u; k c ) must have at least one solution, 
we shall conclude that it has precisely one solution. Let k = k c and suppose there 
exist u x ,u 2 € X, u 1 ^ u 2 , such that P{v}) = h{u l ;k c ) and P(u 2 ) — h(u 2 ;k c ). 
By strict concavity of P we have that P^iu 1 + u 2 )) > \ (P(u 1 ) + P(u 2 )). 
Define u' := ^(u 1 + u 2 ) and note that vl S X. We have that P(u') > h(u'; k c ). 
This implies that there exists k' < k c such that P(u') — h(u';k'). But by 
definition k c is the smallest k for which P(u) = h(u; k) has a solution. We 
arrive at a contradiction and conclude that u 1 = u 2 ; or in other words, that the 
equation P(u) — h(u; k c ) has exactly one solution on X. Denoting this solution 
by u* , it is not hard to see that, at u — u* , the derivative of P with respect to u 
and the derivative of h with respect to u (both of which are defined on the entire 
interval X) must coincide. For suppose f^(u*) < §^(w*), then by continuity 
there exists 8 > such that h(u; k c ) < P(u) for all u such that u — u* < 8. Let 
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u' be one such u. It follows that there exists k' < k c such that P(u') = h(u'; k'). 
This leads to a contradiction and we conclude that h(u; k c ) > P(u). By analogy 
we have that h(u; k c ) < P(u). We conclude that = §£(«*)■ That is, 

1 1 1 N (^Y 

— = — — V (75) 

k c u*N^ I 7T\T 1 ' 



Or equivalently, 



it 1 



-T- 



2V / /„ N 2 , JV 



1^-^)4?:-^ (76) 



3=1 » * ' ^ 1 

Now recall that by definition of u* , we have that 



iV / /r> \ 2 

Equating the right hand side of Eqn. (|76"1) with the right hand side of Eqn. (|77[) 

gives 



This shows that u* is a solution to (|73|) . What remains to be shown is that u* 
is the only solution. Define v, w : X i— » R, 



n, 



and note that, on their respective domains, v is strictly monotonically decreasing 
while w is strictly monotonically increasing. In addition, note that there exist 
a, b G X such that t>(a) > w(a) and v(b) < w(b). Hence, by continuity, there 
must exist a point u' € (a,b) C I such that v(u') = w(u'). Monotonicity of v 
and w implies that this point is unique. It follows that u* is the unique solution 
of J73J) on 1. And by ([77]) we have that 

(79) 



This concludes the proof. 
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Based on the result of Theorem [3] we define the map K : VTR \{0} 

K (0) = U * , (80) 



_L \^ iv h — I Ul 
n 2jj=i Y V 11 * 

where, as before, u* denotes the unique solution of ([73]) on X, given fi. Note 
that, given any realization of u> such that Vto ^ 0, we have that k c — K(Vus). 
We have the following corollary. 

Corollary 7 

1. For all n G M^, f2 ^ 0, we have that < K(O) < 2 1 1 1 1 ^ ; 

,2. i/iere exists f2 £E suc/i that K (Q) = 211171100 if and only if N is even; 

3. for every e > there exist an positive integer N and € VM. N such that 
\K(Q)-\\Q\\ 00 \<e. 

Proof: (Part 1). We show that for all ^ 0, the solution u* of equation (|73[) 
satisfies u* < \/2||0||oo- The result then follows easily. Let u' := v2||^||oo- 
Then we have that 



2±YJl-(*k) 2 > ±Y 1 forall U>U ' 



It follows that 



Km < % < ^ n 



•a 7 ^i=i Y 1 \u> 

for all £1 S KR^. The lower bound if(f2) > ||0||oo was obtained earlier in 
Section [U 

(Part 2.) From the above it follows that K(&Y) — 21101100 if and only if 



2- , ,. - , 
iV ^ V W 7 # • 

j=i ' j=i * 1 — I --f 



or cquivalently, Of = for all It is easy to see that this latter condition is 

never satisfied when N is odd (keeping in mind that J^j % = 0) • Now suppose 
TV is even and pick any c ^ 0. Define 



a, 



c i — 1, 2, ... 2 , 
-c i = 4 + 1,..., N 
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Then we have that € VR N . Moreover, Qf = fi| = c 2 for all It follows 

that K(n) = 2IIOHOO. 

(Part 3.) Let e > be given and suppose N is odd. Pick c ^ and define 



i = l,2,...,JV-l; 
c i = JV. 



Then (f?3")) evaluates to 

2(iV-l) + 2^1-(£) 2 = (iV-l) + " 



(81) 



1 



and it is not hard to see that as N tends to infinity, the solution u* of (fBlj) tends 
to c. Indeed, for N > 2 we have 

(^) 2 = \(N~1) (-\(N-l) + \V(N-ir+8y (82) 

Let ei > 0, and pick TV such that i < e\ and u* < (1 + (\)c. It follows that 



iV 



and hence 
K{Q) := 



1-ei 



(83) 



c. (84) 



Now let ei be given as t\ := ^rpp and choose N accordingly. It follows that 
K(p,) < ||fi||c3o + £■ This concludes the proof. 

We are now ready to present our algorithm, which, given f2, will com- 
pute u* with user-defined precision e > in a finite number of iterations, 

„ =riog3 (JlM-)1 + l. 

Algorithm 1 

1. a := ||fi||oo, b-= \/2||fi||oo 

2. While (b-a)> e, 



3. u:=\{b-a). 



I If 
5. End. 



E,* V I 1- „ ) > \T 



then a := it, eZse := u. 



Once we have an estimate u of it*, we can use (|74j) . replacing u* with m, to 
estimate k c . 
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7 Numerical Example 



We illustrate the results presented in this paper by means of a numerical exam- 
ple. We consider two systems with N — 20 and N = 200 oscillators respectively, 
with frequencies {fif^} as depicted in Figure [5l The frequencies in this example 



1.5 



20 



2.5 

T 

ft 200 





-2.5 



y 



200 



(a) N = 20 



(b) N = 200 



Figure 5: The vector of frequencies Qf := u> 



N e {20,200}, used 



in this example. The natural frequencies u)f were sampled from a normal dis- 
tribution with zero mean and unit variance and relabelled in such a way that 



N 



< i , N 



were sampled from a normal distribution with zero mean and unit variance and 
relabelled such that lo^ < cu^ ... < u>^ (note that this can be done without 
loss of generality). For this particular realization of ui 20 (u> 200 ), we have that 
II^Hoo = 1.7858 (||fi 200 ||oo = 2.3893) and 

^Ef - (Jt) 2 = 0-8015 (0.9139). 

It follows from Corollary © that k c < 2.2281 (2.6145) and by §TZ$, we have 
that k c > = 1.7858 (2.3893). Figure shows the time evolution of the 

magnitude squared of the order parameter, R 2 (t) (previously denoted as L(t)), 
for two different initial conditions and two values of the coupling coefficient, 
k = 2.3 and k = 2.65 (k = 2.1 and k = 2.3). We observe that when k is slightly 
greater than the known lower bound on k c , the value of R 2 (t) converges to a 
constant and inspection shows that the solution x(t) of the system (|T7|) tends to 
a fixed point. On the other hand, when the coupling coefficient is slightly below 
the known upper bound on the critical coupling, the trajectories x(t) appear 
not to converge. Note that in this case we do not know whether the system (fl7|l 
has a fixed point or not, as the condition stated in Corollary [6] is only sufficient 
while at the same time the respective coupling strengths exceed their known 
lower bounds (1.7858 and 2.3893 respectively). To gain more insight into this 
situation let us consider the case N — 20 in some more detail. We fix the 
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Figure 6: Case N — 200: time evolution of the magnitude squared of the order 
parameter, R 2 (t), for two different initial condition (indicated by a dashed and 
solid line respectively), and two values of k. In the left panel, the value of k 
(2.4) is (well) below the known upper bound on k c (2.6145)and the system does 
not converge to a fixed point; in the right panel the value of k (2.65) is slightly 
above the known upper bound on k c and the system converges to a fixed point, 
as expected. 



coupling coefficient at k = 2.1, and numerically evaluate the function P 20 (k, •), 



1 20 

for several values of (3 in the interval [i HQ 20 ^, 1]. We repeat the same com- 
putation for k — 2.3. The result is shown in Figure [51 We observe that the 
equation P 20 (k(3) = (3 does not have a solution on the interval [^||fi 20 ||oo, 1] 
when k = 2.1, but does have a solution when k = 2.3. 

We use Algorithm [T] to compute the 'exact' value of the critical coupling to 
the fifth significant digit. We find that k c = 2.2198 for the case N = 20 and 
k c = 2.6144 for the case N = 200. Note that in both cases, but particularly the 
latter, the upper bounds (2.2281 and 2.6145 respectively) provide good estimates 
of the true values of the critical coupling. 







1 - 







(85) 



8 Conclusion 

We derived necessary and sufficient conditions for the existence of fixed points in 
a finite system of coupled oscillators. In particular, we derived an easy sufficient 
condition in terms of the individual oscillator frequencies (Corollary [3|), which 
we used to compute an upper bound on the critical coupling (Corollary EJ). 
We showed that when no prior knowledge of the distribution of frequencies is 
available, we can still bound the critical coupling in terms of the infinity norm 
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5 10 15 20 5 10 15 20 

t-» t-> 
(a) k = 2.1 (b) k = 2.3 

Figure 7: Case TV = 20: time evolution of the magnitude squared of the order 
parameter, R 2 (t), for two different initial condition (indicated by a dashed and 
solid line respectively), and two values of k. In the left panel, the value of k (2.1) 
is slightly below the known upper bound on k c (2.2281)and the system does not 
converge to a fixed point; in the right panel the value of k (2.3) is slightly above 
the known upper bound on k c and the system converges to a fixed point, as 
expected. 



of the frequencies with their mean removed (Corollary [7]). These bounds were 
shown to be the tightest possible, in the sense that we can find realizations of the 
intrinsic frequencies for which the upper bound is attained, and others for which 
the critical coupling is arbitrarily close to the lower bound. Finally, we proposed 
an efficient algorithm (Algorithm!!} for computing the critical coupling to within 
arbitrary bounds in a finite number of steps. In future work we shall seek to 
extend the present analysis to complex networks of arbitrary topology, and 
investigate more closely the impact of the shape of the distribution of intrinsic 
frequencies on the value of the critical coupling. We shall also consider the 
important question of stability, and present analytical results for the limit case 
when the number of oscillators tends to infinity. 
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